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ABSTRACT 

A supermassive black hole in the nucleus of an elliptical galaxy at the centre of a cool-core 
group or cluster of galaxies is immersed in hot gas. Bondi accretion should occur at a rate 
determined by the properties of the gas at the Bondi radius and the mass of the black hole. 
X-ray observations of massive nearby elliptical galaxies, including M87 in the Virgo clus- 
ter, indicate a Bondi accretion rate Mb which roughly matches the total kinetic power of the 
jets, suggesting that there is a tight coupling between the jet power and the mass accretion 
rate. While the Bondi model considers non-rotating gas, it is likely that the external gas has 
some angular momentum, which previous studies have shown could decrease the accretion 
rate drastically. We investigate here the possibility that viscosity acts at all radii to transport 
angular momentum outward so that the accretion inflow proceeds rapidly and steadily. The 
situation corresponds to a giant Advection Dominated Accretion Flow (ADAF) which extends 
from beyond the Bondi radius down to the black hole. We find solutions of the ADAF equa- 
tions in which the gas accretes at just a factor of a few less than Mg. These solutions assume 
that the atmosphere beyond the Bondi radius rotates with a sub-Keplerian velocity and that the 
viscosity parameter is large, a > 0.1, both of which are reasonable for the problem at hand. 
The infall time of the ADAF solutions is no more than a few times the free-fall time. Thus 
the accretion rate at the black hole is closely coupled to the surrounding gas, enabling tight 
feedback to occur We show that jet powers of a few per cent of Mbc^ are expected if either a 
fraction of the accretion power is channeled into the jet or the black hole spin energy is tapped 
by a strong magnetic field pressed against the black hole by the pressure of the accretion flow. 
We discuss the Bemouilli parameter of the flow, the role of convection, and the possibility that 
these as well as MHD effects may invalidate the model. If the latter comes to pass, it would 
imply that the rough agreement between observed jet powers and the Bondi accretion rate is 
a coincidence and jet power is determined by factors other than the mass accretion rate. 

Key words: X-rays: galaxies — galaxies: clusters — intergalactic medium: accretion, accre- 
tion discs — black hole physics 



1 INTRODUCTION 

The nuclei of massive elliptical galaxies at the centres of cool-core 
groups and clusters of galaxies have powerful relativistic jets which 
inject energy into the suiTounding hot gas. This prevents the intra- 
cluster gas from radiatively cooling and collapsing onto the galaxy, 
thus stifling its growth (McNamara & Nulsen 2007 and references 
therein). Such feedback is now a common ingredient in our under- 
standing of the evolution of massive galaxies (Croton et al. 2006; 
Hopkins et al. 2006). 

The mode of fuelling of the massive black hole at the galaxy 
nucleus, which energises the jets, is unclear. Bondi (1952) accre- 
tion is often invoked since the black hole is sitting in the dens- 



* E-mail:marayan@cfa.harvard.edu 
t E-mail:acf@ast. cam. ac.uk 



est part of the hot cluster (or group) atmosphere. Observations of 
the gas around the Bondi radius in M87 indicate that Bondi ac- 
cretion may indeed provide a suitable mass supply rate (di Mat- 
teo et al. 2003). Others argue that it cannot provide enough fuel 
to power more powerful, distant objects (Rafferty et al. 2006), and 
that cold gas clouds may instead be required (Pizzolato & Soker 
2005). Regardless, the Bondi model considers gas with vanishing 
angular momentum, whereas in a realistic situation the incoming 
gas is likely to have non-negligible rotation. Hence, it is not clear 
that the Bondi accretion rate Mb is at all relevant. 

Most nuclei in the centres of cool core clusters show no sign 
of a dense, radiatively-efficient accretion disc. Some of the most 
powerful ones do not even show any detectable X-ray point source 
(Hlavecek-Larrondo & Fabian 2011), which is difficult to explain 
in cold mode accretion. In the case of M87, there is clear evi- 
dence that both the accretion flow and the jets themselves are ra- 
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diatively inefficient (di Matteo et al. 2003). This indicates that the 
flow must be advection dominated, i.e., the gravitational energy re- 
leased in the flow must be carried into the centre rather than ra- 
diated locally (see Narayan, Mahadevan & Quataert 1998; Kato, 
Fukue & Mineshige 2008; Narayan & McClintock 2008; for re- 
views of advection-dominated accretion flows, or ADAFs). A frac- 
tion of the energy must then be efficiently transfeiTed to the jets 
once the accreting gas reaches the centre. 

The above conclusion is supported by a study of 8 other mas- 
sive nearby elliptical galaxies where the gas properties close to the 
Bondi radius can be observed or reasonably extrapolated (Allen et 
al. 2006). In all these cases, the Bondi mass accretion rate Mg deter- 
mined at the Bondi radius tb correlates well with the power of the 
jets P^, where the latter is measured from the bubbles inflated by the 
jets in the surrounding gas. Writing the jet power as = rj^M^c^, 
the jet production efficiency factor 77j is found to be about 2 per 
cent. This is a rather large efficiency and underscores the need for 
mass at approximately the Bondi accretion rate reaching the gravi- 
tational radius of the black hole rg. There is little room for any in- 
efficiency in the transport of mass to the centre, e.g., through mass 
loss in outflows along the way. 

We are concerned here whether an ADAF can be established 
in galactic nuclei and whether the mass accretion rate is comparable 
to the Bondi rate. The range of jet power in the systems discussed 
above is between 10"*^ — 10*^ ergs"' , so the Eddington ratio (power 
emitted to Eddington limit) is 10~^ — 10^^ for a black hole of mass 
lO^Mcv, and ten times less for lO'^M^. This is very much in the 
regime where an ADAF is expected (Narayan & Yi 1995b; Narayan 
& McClintock 2008). Moreover, as noted by Narayan & Yi (1994, 
1995a) and Fabian & Rees (1995), and confirmed in more recent in- 
vestigations (Narayan & McClintock 2008), the large thermal pres- 
sure of an ADAF may be especially good for the production and 
coUimation of jets. Thus, it is natural to consider an ADAF-like 
accretion model for systems with powerful jets. 

ADAFs have been well studied since the work of Narayan 
& Yi (1994, 1995ab) and Abramowicz et al. (1995). However, in 
much of the previous work, the outer edge of the solution was 
generally taken to be either of a self-similar form (e.g., Chen, 
Abramowicz & Lasota 1997; Popham & Gammie 1998) or a ge- 
ometrically thin disk that evaporates to form the ADAF (e.g., 
Narayan, Kato & Honma 1997; Manmoto et al. 2000). Neither of 
these boundary conditions is relevant for understanding accretion 
from an external medium. Since an ADAF is essentially space fill- 
ing, we expect the accretion flow to match smoothly on to the ex- 
ternal medium without any shocks or other kinds of discontinuities. 
We investigate in this paper exactly how this matching occurs when 
we have a slowly rotating external medium (eq. |23]gives a quanti- 
tative measure of what we mean by slow rotation). 

Previous studies of Bondi-like accretion with angular momen- 
tum have generally considered inviscid flows. Proga & Begelman 
(2003) carried out two-dimensional axisymmetric simulations and 
showed that an equatorial torus forms because of the angular mo- 
mentum barrier and that this torus constrains the amount of polar 
accretion. Krumholz, McKee & Klein (2005) extended their work 
and developed approximate formulae for the mass accretion rate 
as a function of the vorticity of the external gas, and Cuadra et al. 
(2006) carried out detailed simulations of inviscid accretion on to 
Sagittarius A* at the Galactic Centre. Recently, Inogamov & Sun- 
yaev (2010) proposed an accretion model for M87. As in the other 
studies cited here, the centrifugal baiTier causes the inviscid ac- 
creting gas to form a torus well inside rg. Inogamov & Sunyaev 
assume that viscosity then turns on at smaller radii and suggest that 



the torus will thus feed a standard thin accretion disc on the in- 
side, which might evaporate into an ADAF at yet smaller radii. The 
presence of the thin disk segment in their model causes the total 
inflow time of the gas from the Bondi radius to the black hole 
gravitational radius rg to be far longer than for a Bondi flow or (as 
we shall see) an ADAF. Self-adjustment of the feedback, in which 
the jet power responds to conditions (e.g. cooling time) beyond r^, 
then becomes very difficult, with large hysteresis expected. 

In contrast to the above studies, we are interested in viscous 
accretion. The closest paper to our work is Park (2009). For tech- 
nical reasons, that work focused on extremely hot external media 
(Text > 10^ K) for which the Bondi radius is much closer to the 
black hole than in real systems. We consider more realistic exter- 
nal conditions (Text ~ 10^^^ K). We also study in more detail the 
transition from a Bondi flow to an ADAF as the external rotation is 
varied. 

As in Park (2009), we require the flow to be continuous out 
through and beyond. Such a model ensures that the accretion 
power is as well coupled with the conditions in the outer gas 
as possible, thereby allowing for the most efficient feedback. We 
moreover require that outflows, and significant radial exchanges of 
energy within the ADAF, are suppressed. We postulate that rela- 
tivistic jets are created and mechanically powered very efficiently 
(but very radiatively inefficiently) by the accreting gas close to the 
black hole, but how this occurs is beyond the scope of the present 
work. We limit ourselves to a more basic question: Can an idealised 
ADAF transfer a high enough mass accretion rate from beyond 
down to rg? 



2 SPHERICAL ADAF MODEL 

2.1 Viscous Accretion Flow: Conservation Laws 

Since we are primarily interested in slowly-rotating, steady, viscous 
accretion flows, we assume that the density and pressure of the gas 
are distributed spherically at each radius. We also assume that all 
quantities are independent of time (steady state assumption). We 
thus focus only on radial variations. Under these assumptions, the 
mass accretion rate M at radius r is given by 



M = —4nrpv — constant, 



(1) 



where p{r) is the density and v(r) is the radial velocity; the latter 
is taken to be negative when gas flows inwards. When consider- 
ing accretion flows in which rotational support is important, e.g., 
geometrically thin disks, or ADAFs with more rotation than we 
consider here, the factor 4jzr^ in the above relation is replaced by 
{2Kr)(2H), where H{r) is the "vertical" scale height of the gas at 
radius r. In the simpler approximation considered here, we effec- 
tively set H = r, which could be interpreted as a geometrically very 
thick disk. Except for this difference, the equations we consider are 
identical to those described in Narayan et al. (1997). 

To mock up relativistic gravity in our Newtonian model, we 
assume a gravitational potential (Paczynski & Wiita 1980) 



GM 



2GM 



(2) 



where M is the mass of the central black hole. Correspondingly, the 
Keplerian angular frequency Q.]^ is given by 



0.i = 



GM 



{r-rgYr 



(3) 
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Making the replacement p = pc^, where Cj is the (isothermal) 
sound speed, we write the steady state radial momentum equation 
as 



(4) 



where fi is the angular velocity of the gas on the equatorial plane. 
Our spherical model is most accurate when the centrifugal acceler- 
ation on the gas is much weaker than the gravitational acceleration; 
this corresponds to the condition Q? <C fi|.. 

We model viscosity via the standard a-prescription (Shakura 
& Sunyaev 1973) in which the kinematic coefficient of viscosity v 
is written as 



V = acsH = acsK 



(5) 



with a taken to be a constant. However, we do not set the shear 
stress equal to ap, but use a more physical prescription in which 
the stress is proportional to the angular velocity gradient: 



shear stress = c,.,* = vprdQ./dr. 



(6) 



The angular momentum equation then takes the form (Narayan et 
al. 1997) 



Id/ 5rff2 
pr^ dr \ ^ dr 



which on integration gives 

dSl _ v{£lr^ - j) 
dr ar^c\ 



(7) 



(8) 



The quantity j is an integration constant with dimensions of spe- 
cific angular momentum. 

Finally, energy conservation gives 

pv dc^ 2 dp ^ f dQ.\^ 

where y is the adiabatic index of the gas, which is set to 5/3 for all 
the numerical models presented here. The left-hand side of equa- 
tion l|9} represents the Lagrangian time derivative of the entropy 
of the gas. This term is usually referred to as the energy advec- 
tion term. The term on the right-hand side of the equation describes 
the heating rate due to viscous dissipation. In the spirit of a radia- 
tively inefficient flow (ADAF), we ignore radiative cooling alto- 
gether, Thus, we set advection equal to heating to obtain equation 

We should note the following inconsistency in the above equa- 
tionfl While we have included the effect of viscosity through the 
shear stress in the angular momentum equation (|7), we have ne- 
glected corresponding terms in the radial momentum equation 
Under the assumptions of our model (pure radial flow, no gradi- 
ents in the transverse direction, etc.), the rr component of the stress 
takes the form (Landau & Lifshitz 1959) 



2 4 d(v,/r) , 
Grr = -Pci + -npr '^'^ + ^p 



1 d{r\'r) 



(10) 



where ^ is the kinematic bulk viscosity. The last term in equation 
(HI should thus be written as ( 1 /p)d{arr) /dr with the above form 
of <7rr, not just as —{i/p)d{pc1)/dr. 

Traditionally, in accretion disk models, the viscous terms in 
Grr are neglected on the grounds that v,. is much smaller than c, 
and so these terms are small compared to the pressure. This is 

' We thank the referee for pointing this out to us 



no longer obvious for the slowly-rotating solutions presented here, 
for which Vr is fairly large. Nevertheless, we make this assump- 
tion for easy comparison with previous work. A major goal of the 
present work is to study the transition from the rapidly-rotating 
ADAF regime to the non-rotating Bondi regime. The viscous terms 
in Grr survive even for pure radial flow and ought to be included 
in a self-consistent model of spherical inflow. Since these terms are 
neglected in the Bondi model, in the same spirit we neglect them in 
our model as well. We leave for the future an investigation of the 
effect of these terms on both the Bondi solution and our slowly- 
rotating solution. 



2.2 The Inner Supersonic Region 



The equations in i|2. 11 correspond to a viscous rotating accretion 
flow. Once the accreting gas passes inside the sonic radius and 
becomes supersonic, we expect viscosity to be much reduced and 
perhaps even to vanish (Narayan 1992; Kato & Inagaki 1994; Kley 
& Papaloizou 1997). For this region of the flow, we simplify the 
equations by setting a = 0, thus dropping all terms related to vis- 
cosity. From the angular momentum equation (|7} we see that the 
specific angular momentum is then a constant: 

r < : = Q.r^ — constant. (11) 

Similary, from the energy equation (|9), we see that the entropy of 
the gas is constant: 



r < : s-,„ = 



,(7-1) 



: constant 



(12) 



Finally, by combining the various conservation laws, we can show 
that the Bernoulli parameter ^ of the gas is also constant. This 
gives the condition 



r<rs: 



v2 el GM 



2 2f-2 {r — rg) 

(y- l)r2(r-l) Ivllr-l) \^4;;: J 

constant. 



M \ 



(7-1) 



(13) 



Using the final relation, along with the values of the conserved 
quantities M, jjn and .^i", we can solve for the radial velocity v 
as a function of r in the supersonic region. This immediately gives 
all the other quantities. 



2.3 Boundary Conditions 

Our model accretion flow consists of two regions: a viscous sub- 
sonic region which extends from the sonic radius out to some 
large outer radius rout, and an inviscid supersonic region which ex- 
tends from the sonic radius down to the black hole. Finding the 
solution in the viscous region requires solving a boundary value 
problem involving a number of differential equations^ Equations 
(0, {8} and ([9} represent three first order ordinary differential equa- 
tions, which require three boundary conditions. In addition, the 
constants M and j are eigenvalues, which require two more bound- 
ary conditions. Finally, the location of the sonic radius has to 



^ Once we have the solution in the viscous subsonic region, we can com- 
pute the values of M, ii„ and ^ at r,. We can then directly calculate the 
solution in the supersonic region. The latter involves only algebraic equa- 
tions (see i|2.2t , not differential equations. 
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be determined as part of the solution, so this requires yet another 
boundary condition. Thus we need to supply a total of six boundary 
conditions. 

The three differential equations (|4j, (|8) and (|9), in combina- 
tion with equation can be reduced to the following relation: 



(yc 



,2 2^^1n^i 



dr 



2ycl , {y-\){ar^-j)\ 



,(14) 



which becomes singular when yc^ — = 0. The radius at which 
this happens is the sonic radius r^, where the flow speed |v| is equal 
to the adiabatic sound speed Cg^/y. In order to have a smooth flow 
through rs, the quantity on the right hand side of equation J14t 
should vanish. We thus obtain the following two boundary con- 
ditions: 

r = rs : yc^-v^ = 0, (15) 

r = rs ■ (Llj( — Ll)rs 1 5 = U(lo) 

I's Olr^Cs 

Viscous accretion flows have another boundary condition on 
the inside, which is usually applied as a no-torque condition at 
some radius In the most elaborate version of the theory, one would 
apply the no-torque condition at the black hole horizon {r = 
however, this tends to make the numerical computations very dif- 
ficult. It also introduces some subtlety into the problem since the 
behavior of viscosity in the supersonic plunging region of the flow 
is poorly understood (Narayan 1992; Kato & Inagaki 1994). In this 
paper, we have assumed for simplicity that viscosity vanishes inside 
the sonic radius. One consequence of this approximation is that the 
specific angular momentum of the gas £lr^ becomes constant in the 
plunging region. Motivated by this fact, we set d{0.r^)ldr = as a 
boundary condition on the viscous solution at r = rj. This ensures 
a smooth transition across the sonic radius. Making use of equation 
([8j, the condition can be written as 



(17) 



The remaining three boundary conditions are applied at the 
outer radius r^ui of the solution (Yuan 1999). We choose rout to be 
large enough that it lies well into the external uniform medium. In 
analogy with the Bondi problem, the temperature of the external 
gas, or equivalently the sound speed, and the density of the gas 
provide two outer boundary conditions: 

^ = '"out ■ Cj = Cout, (18) 

'■ = ''out : P=Pout- (19) 

In the numerical solutions presented here, we set Cout = lO^-'c, 
which corresponds to a temperature of 6.5 x 10* K (assuming a 
mean molecular weight of 0.6), a reasonable choice for the inter- 
stellar medium at the centre of a galaxy. In the case of the density, 
we arbitrarily select pout = 1- After the fact, we can rescale the den- 
sity profile by a constant factor so as to satisfy the required value 
of Pout- This approach is allowed by the fact that the equations are 
linear in the density^ 

^ This boundary condition is needed only for the more physical viscous 
stress prescription (eq. |6} used here. If the shear stress is written in the 
simpler form a p, there is one fewer differential equation and the additional 
boundary condition is not needed (see Narayan et al. 1997 for a discussion). 
In fact, since pressure never vanishes in an accretion solution, the a p stress 
prescription does not have vanishing stress at any radius. 

This is true only because we have ignored all cooling terms. If we include 
a detailed cooling model, the physics will no longer be linear in p . 



For the third boundary condition, we fix the angular velocity 
of the external gas: 



r = rout : fl = flo 



(20) 



However, we note the following complication. Because we are 
solving viscous accretion equations with a constant a, the solu- 
tion naturally tends to a state of rigid rotation on the outside. For 
radii outside the Bondi radius. 



GM 



'^"out 



(21) 



the black hole gravity is too weak to influence the dynamics of 
the gas - pressure is more important here. As a result, viscosity 
drives the gas towards dO-jdr = 0. In a real galactic nucleus, this is 
precisely the region where the gravitational potential of the galaxy 
will take over and the gas will transition to the rotation curve of the 
galaxy (see Quataert & Narayan 2000 for a study of Bondi accre- 
tion in such a potential). Since we have not included the galactic 
contribution to the potential (|2j, our equations enforce a rigidly 
rotating external medium. The problem with this is that the cen- 
trifugal acceleration fl^r increases without bound at large radius, 
which is unphysical. To avoid this problem we choose rout to be 
only a factor of a few (not more than 10) larger than r^. By making 
this choice, we ensure that the centrifugal acceleration does not be- 
come too large on the outside. At the same time, we make sure that 
rout is large enough for the solution to asymptote to the conditions 
in the external medium. 

The parameter iiout determines whether the external gas is 
rotating slowly or rapidly. The boundary between the black hole- 
dominated accretion flow and the external medium is located at 
r ~ rs, and the Keplerian angular velocity fl/j-^s at this radius is 
given by 



■a:,b 



GM 



1/2 



= 2 



(^'out 



-3 



(22) 



We thus define the following dimensionless rotation parameter iM: 
^^;^ = liloutf — V- (23) 

^l-KS) ^ \Cout/ 

When ^ <C 1, we say that the external medium rotates slowly, 
whereas as ^ approaches unity, the medium rotates rapidly. We 
are primarily interested in the slowly rotating case. 

If the external medium rotates slowly enough, the gas may 
be able to accrete directly into the black hole without any need 
for viscous transport of angular momentum. We would then have 
something very similar to the Bondi solution. The critical angular 
momentum of the external gas at which we expect this transition 
to take place is the specific angular momentum of the marginally 
stable orbit f^s, which for the potential (|2j is 



(24) 



Correspondingly, we can express the angular momentum of the ex- 
ternal gas as the following dimensionless ratio 



^^^ = ^=0.I36I^outf— 



(25) 



When ^ ^ 1, we expect the flow to be viscously driven and to 
resemble an ADAF solution, whereas when ^ <C 1, the flow should 
be practically identical to the Bondi solution. These expectations 
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Figure 2. Left: Shows the location of the sonic radius rj as a function of 
the angular momentum parameter for solutions with a = 0.1, 7= 5/3, 
Cout = lO^^c. The upper dotted line indicates the sonic radius for a pure 
non-rotating Bondi solution, and the lower dotted line shows the radius of 
the marginally stable orbit r„^s - Note the sudden transition from a Bondi-like 
flow for ^ < 1 to a rotation-supported ADAF for >2. Right: Shows the 
corresponding mass accretion rates M in units of the Bondi accretion rate 
Mb . The mass accretion rate is only a factor of three smaller than the Bondi 
rate even when ^ is as large as ~ 10^. 



Figure 1. Representative solutions of the model equations for a = 0. 1, 
7= 5/3, Cout = lO^^c. The four solutions shown have (flout, ''j) = 

(0.624 X 10"', 0.31, 85, 3.436); (0.851 x 10"'", 0.043, 12, 3.663); 
(0.132 X 10"'°, 0.0066, 1.8, 142.0); (0.831 x 10"'^, 0.00042, O.ll, 416.7), 
respectively. The solid dots indicate the positions of the sonic radii and are 
helpful for identifying the solutions. In addition, a few curves are labeled 
by their values of The vertical dotted Hnes correspond to the location of 
the Bondi radius (eg. 1211 . and the sloping dotted line in the top left panel 
shows the Keplerian angular frequency Q.k . 

are borne out by the numerical solutions described in ij3] For our 
choice of tout = W^^c, Jff = I corresponds to = 0.0037. 



3 NUMERICAL RESULTS 

Since the viscous accretion equations tend to be very stiff, we 
use a relaxation method (Press et al. 1992) to solve them0 Fig- 
ure 1 shows sample solutions corresponding to a = 0.1, y = 5/3, 
Cout = 10"^ and po^t = 1 (the value of Pout is arbitrary since we can 
rescale the density profile to any external density as needed, i]2.3b . 
Four solutions are shown, corresponding to .if = 85, 12, 1.8, 0.11, 
respectively (compare with Fig. 1 in Park 2009). Note that the rota- 
tion parameter 1% is small for all the solutions, so these truly repre- 
sent slowly-rotating flows. Even the most rapidly rotating solution 
(,^ = 0.31) has a centrifugal support of only 10% of Keplerian at 
r = rB. 

The solution with ^ = 0.11 - the lowest curve in the top- 
left panel of Fig. 1 - is clearly in the Bondi regime since the gas 
has negligible outer specific angular momentum relative to /ms- The 
sonic radius r^, shown by the black dot, is located at 417 rg, which is 
almost exactly where a pure non-rotating Bondi flow has its sonic 
radius for our choice of <S(r), Cout and y. The two solutions with 
^ = 85 and 12 (the highest two curves) are definitely rotation- 
dominated. The gas in these solutions has too much angular mo- 

^ The simpler shooting method is adequate if the outer radius is not 
too large, e.g., ra^/rg < 10^. However, for realistic external media with 
Cout/c ~ 10"^, we need to calculate solutions out to Tom/'g > lO*". In our 
experience, relaxation is the only sure way to obtain such solutions. 



mentum to permit steady accretion in the absence of viscosity, so 
the accretion flow settles down to a viscously-driven ADAF solu- 
tion. Correspondingly, the sonic radius is close to the marginally 
stable orbit, r^s = 3rg. The solution with ^ = 1.8 represents a 
transition state between the Bondi and ADAF regimes. Its sonic 
radius is at an intermediate location, = I42r^. 

In Figure 2, the top-left panel shows how the sonic radius 
moves as we change Jf. For all values of ^ < 1, is located at 
the position one would calculate for the non-rotating Bondi prob- 
lem (upper dotted line), while for Jf greater than a few, is close 
to (lower dotted line). The transition between these two regimes 
is quite sudden, with most of the change happening over the range 
1.5<^<2. 

The bottom two panels in Fig. 1 show the profiles of density p 
and pressure p = pel for the same four solutions as in the top left 
panel. Even though the rotation profiles of these solutions are very 
different, and their sonic radii move around considerably, the pro- 
files of p and p are nearly identical. The insensitivity to the location 
of is at least in part because we selected y = 5 /3 , which is known 
to be a critical value of the adiabatic index both for the Bondi prob- 
lem and for ADAFs. Nevertheless, it is clear that in many respects, 
an ADAF is very similar to a Bondi flow. 

The top-right panel in Fig. 1 shows the radial velocity pro- 
files of the four solutions. We see that the radial velocity is smaller 
for a rotating ADAF (the solutions with ,5f = 85, 12) compared 
to a slowly-rotating Bondi-like flow {J^ = 0. 1 1). Since the density 
profiles of both kinds of solution are nearly the same, this means 
that the mass accretion rates are different. This is illustrated in the 
panel on the right in Fig. 2 (compare with Fig. 2 in Park 2009), 
which shows that M decreases as the rotation of the external gas 
increases. The effect is quite modest, however — the total range of 
M in our solutions is only a factor of three, though this is an artifact 
of our choice of a relatively large value of a = 0. 1. 

To illustrate the effect of a. Figures 3 and 4 show results cor- 
responding to ^ = 13.5 (.^ = 0.05) and six values of a: 0.316, 
0.1, 0.0316, 0.01, 0.00316, 0.001. The rotation profiles are nearly 
the same for different values of a, with only small variations. 
More interesting is the behavior of the sonic radius, which moves 
in a very systematic way as a is varied. For a = 0.316, we find 
rs = 17.5, which is well outside the radius of the marginally sta- 
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Figure 3. Solutions for 7=5/3, c^m = W c, ilout = 10"'", = 0.05, 
if = 13.5, and six values of a: 0.316, 0.1, 0.0316, 0.01, 0.00316, 0.001. 
The six solutions have r, = 17.54, 3.656, 2.694, 2.335, 2.174, 2.114, re- 
spectively. The solid dots indicate the positions of the sonic radii and are 
helpful for identifying the solutions. In addition, a few curves are labeled 
by their values of a. The vertical dotted lines correspond to the location of 
the Bondi radius , and the sloping dotted line in the top left panel shows 
the Keplerian angular frequency 01^- 



ble orbit = Sr^. With decreasing a, rs moves in until it is well 
inside rms. The pattern is very similar to that seen in the ADAF 
models described in Narayan et al. (1997) and noted in many other 
papers (e.g., discussion of the slim disk model by Abramowicz et 
al. 2010). 

The radial velocities shown in Fig. 3 decrease proportional to 
a, and so do the mass accretion rates (Fig. 4 right panel). This is 
consistent with the predictions of the analytical ADAF model de- 
scribed in Narayan & Yi (1994), and is in qualitative agreement 
with Park (2009). In particular, we agree with Park's conclusion 
that low angular momentum flows resemble Bondi accretion, and 
that their accretion rates approach the Bondi rate Mb (which is 



equal to 3.65 x 10 in our units where rg 



1) as the ex- 



ternal rotation decreases. However, we also see some quantitative 
differences. Park suggests on the basis of his numerical solutions 
that the accretion rate scales approximately as M/Mq ~ 9a/.^ (our 
parameter is the same as X in Park's notation). We do not repro- 
duce the scaling with J^. For instance, Fig. 2 shows that, at fixed 
a = 0. 1, M changes by only a factor of '-^ 3 as ,^ changes by nearly 
a factor of 100. One reason for this difference could be that we have 
considered solutions with r^/rg = 10^'^, whereas Park's solutions 
are closer to 10^. A more thorough exploration of solutions in the 
three-dimensional space of a-M-{rB / rg) would be worthwhile to 
map out how the accretion rate varies with these parameters. 

The bottom two panels of Fig. 3 show that, with varying a, the 
density and pressure are largely independent of a, just as we earlier 
found them to be independent of flout - The insensitivity of the cen- 
tral pressure to any parameter other than the external density Pout 
and sound speed Cout is a strong result and may have implications 
for jets 



Figure 4. Left: Shows the location of the sonic radius as a function of 
the viscosity parameter a for the solutions described in Fig. 3. The dotted 
line indicates the radius of the marginally stable orbit r^s- Large values of 
a cause the sonic radius to move outside r^s while small values of a have 
the opposite effect. Right: Shows the corresponding mass accretion rates M 
in units of the Bondi accretion rate Mb . Note that M is a steep function of 
a, varying almost linearly. This is expected for an ADAF. 



4 DISCUSSION 

The primary goal of this paper was to estimate the rate at which 
mass accretes on to a supermassive black hole from rotating gas in 
the nucleus of a galaxy. We find that the answer depends on both 
the dimensionless rotation parameter of the gas ^ (ea.l23t and the 
viscosity parameter a (eq.[5}- For fixed a = 0. 1, the accretion rate 
is within a factor of a few of the Bondi rate for all values of (Fig. 
2), i.e., for this value of a accretion is nearly as efficient in the 
presence of rotation as in its absence (the classic Bondi problem). 
Figure 4 shows the variation of M for fixed = 0.05 (a reasonable 
value, see Inogamov & Sunyaev 2010) and different values of a. 
Here the variation is much larger. The accretion rate is suppressed 
by a large factor when a <C 1. Hence, the answer to our primary 
question depends very much on the value of a. 

King, Pringle & Livio (2007) have examined a variety of ob- 
servational evidence and conclude that many observed accretion 
disks have a --^ 0. 1 — 0.4. In addition, Sharma et al. (2006) argue 
that ADAFs have, if at all, even larger values of a compared to 
standard thin disks. Thus, the models with a = 0.316 and 0.1 in 
Fig. 4 may be the best match to real radiatively-inefficient accre- 
tion flows in galactic nuclei. If so, we can expect viscous accre- 
tion via an ADAF to be quite efficient in galactic nuclei: Madaf ~ 
(0.3- 1) X Mb- 

As discussed in ^ our picture is that a fraction Tja^c of the 
accretion energy near the black hole is somehow converted to jet 
mechanical energy. If the accretion rate is equal to the Bondi rate, 
then observations require about 2 per cent of Mbc^ to be trans- 
ferred to the jet. From the estimates given above, we see that accre- 
tion of rotating gas via an ADAF requires an efficiency of perhaps 
'7acc 5%. While we do not have a model of how the jet is actually 
launched, an efficiency of 5% does not appear implausible. 

It should be emphasized that the present study differs from 
most previous discussions of this problem in the literature in that 
we consider a steady viscous flow extending from beyond the 
Bondi radius down to the black hole horizon. Viscosity enables 
our solutions to overcome the centrifugal barrier and to accrete 
steadily, just as in the standard thin accretion disk model (Shakura 
& Sunyaev 1973; Frank, King & Raine 2002). In contrast, most 
other studies of quasi-spherical accretion with rotation (except Park 
2009, see ijlj have considered inviscid accretion. Those results de- 
pend critically on the assumption that accretion is arrested once the 
gas hits the centrifugal barrier. 



© 0000 RAS, MNRAS 000, 000-000 



A Spin on Bondi Flow 7 



From our numerical solutions, we can calculate the time re- 
quired for gas to travel from the Bondi radius down to the black 
hole. For the two models with a = 0.316 and 0.1 in Fig. 4, the ac- 
cretion time ?ADAF for the ADAF solution is no more than twice 
as long as the accretion time in the non-rotating Bondi solu- 
tion. Even for the rapidly rotating solution with ^ = 0.31 in Fig. 
2, fADAF is only ~ 3fBLl This is very encouraging. For a turbu- 
lent external medium, we expect the rotation of the external gas 
near the Bondi radius to vary on a time scale of the order of tens 
of fB (assuming a turbulent Mach number --^ 0. 1, Inogamov & Sun- 
yaev 2010). Since our ADAF solutions have an accretion time much 
shorter than the turbulence time, there is no problem setting up the 
steady state conditions we assume. More importantly, the short ac- 
cretion time guarantees that any feedback from the ADAF via jets 
will occur rapidly compared to the dynamical time of the exter- 
nal medium. Such instantaneous feedback is generally assumed in 
most current models of feedback. 

While we have focused so far on the accretion flow as the 
source of jet power, a popular alternative hypothesis involves the 
black hole. Blandford & Znajek (1977) developed a scenario in 
which the rotational energy of a spinning black hole is tapped by a 
magnetic field and carried away in a magnetized jet. For a magnetic 
flux <J>B threading the horizon, Tchekhovskoy, Narayan & McKin- 
ney (2010) give a fairly accurate estimate of the jet power in ge- 
ometrized units (GM = c = 1), 

Pj = k^lal, O.H = ^, rH=2M[l + (l-fl,)2]i/2, (26) 

where k ~ 0.05 is a dimensionless number that depends weakly on 
the field geometry, a* = a/M is the dimensionless spin parameter 
of the black hole, and rn and Slf] are the radius and angular fre- 
quency of the black hole horizon. The magnetic flux is given by 
<I>B = ATtrjj\B\H where \B\h is the field strength at the horizon. 
Thus 

Pj = 32;r^fcpinag'-H"* ~ 50pmagrjjal, (27) 

where pmag = |S|^/8;r is the magnetic pressure at the horizon. 

To maintain a magnetic field on the horizon, it is neces- 
sary to keep the field lines in place by means of an external 
pressure. We assume that this pressure is supplied by the accre- 
tion flow. Consider first the Bondi non-rotating solution. In our 
units (rg = c = Pout = 1), the thermal pressure of the Bondi so- 
lution at r = 2rg ij3 Ptherm(2''^) = 3.34 x 10^, and the ram pres- 
sure is pmm{2rg) = pv^ = 6.51 X 10*, giving a total pressure of 
Ptot(2fg) = 6.85 X 10*. We assume that the total pressure is what 
confines the central field and write the magnetic pressure at the 
base of the jet as Pb = i]BPtot(2r(,), where T]^ < 1 is a proportion- 
ality constant. We also replace r/j by rg in equation l l27b . which is 
an overestimate for a rapidly spinning hole but is probably a rea- 
sonable simplification. We then find 

Bondi : Pj 3.4 x 10^ risal t^boIMbc^- (28) 

This interesting result shows that the jet power scales directly with 
the Bondi accretion energy rate Mbc^ and varies strongly with the 
spin of the black hole. Our guess is that tjb is probably in the 
range 0.1 — 1. Therefore, provided the black hole does not rotate 

* For very small values of a, we do we find Iadaf S> tB,e.g., for a = 0.001, 
^ = 0.05, we obtain fADAp/fB = 156. However, such low values of a seem 
unlikely. 

^ We cannot use r = rg because our Newtonian potential l|2) is singular 
there. 




log(r/rg) log(r/r|.) 

Figure 5. Variation with radius of the Bernoulli parameter (eq. ll3K nor- 
malized by the gravitational potential, for the four solutions shown in Fig. 
1 (left panel) and the six solutions shown in Fig. 3 (right panel). The solid 
dots indicate the positions of the sonic radii. Curves are labeled by their 
values of ^ in the left panel and a in the right panel. 

too slowly, a Bondi flow could easily support the jets seen in obser- 
vations. 

We now compute the central pressures in the ADAF solutions. 
The four solutions shown in Fig. 2 have pressures ranging from 
Ptot(2r^) = 1.49 X 10* for = 0.31 to ptot(2r,,) = 6.84 x 10* for 
= 0.00046, while the a = 0.316 and 0.1 solutions in Fig. 4 have 
Ptot{2rg) = 6.27 X 10* and 2.59 x 10*, respectively. These are the 
solutions most relevant for our problem, and their pressures are be- 
tween 20% and 100% of the Bondi pressure. Thus, we obtain the 
following estimate for the jet power in the presence of an ADAF, 

ADAF : Pj ^ (0.2 - 1) x rjBalMBC^- (29) 

The jet efficiency factor in this model is Tjj = (0.2 — 1 ) x rjB and a 
net efficiency of 2% seems quite plausible. 

To summarize, in terms of energetics at least, we have two 
viable mechanisms to power jets via accretion in galactic nuclei: (i) 
by tapping a fraction of the accretion energy, and (ii) by confining 
a strong magnetic field around a spinning black hole and extracting 
energy from the hole. Neither mechanism requires us to postulate 
extreme conditions or to stretch parameters. However, all of our 
results are based on the simple one-dimensional model described 
in Sj2] Unfortunately, there are several important caveats that need 
to be discussed. 

Narayan & Yi (1994, 1995a) showed that the accreting gas in 
an ADAF has a positive Bernoulli parameter ^ (see eq.[T3]for the 
definition). This means that the gas is gravitationally unbound, and 
so these authors suggested that ADAFs would have strong outflows 
and jets. The mass conservation equation ([TJ explicitly ignores such 
outflows. How strong the outflows are is difficult to estimate from 
first principles, though Blandford & Begelman (1999) suggested 
that the effect may be so strong that Madaf at the black hole might 
be reduced by orders of magnitude. 

Numerical hydrodynamic simulations have confirmed that the 
mass accretion rate is indeed suppressed (Stone, Pringle & Begel- 
man 1999; Igumenshchev & Abramowicz 1999, 2000; Igumen- 
shchev 2000), though there is no consensus on whether gas truly 
escapes to infinity, or if the Bernoulli parameter is even relevant 
(Abramowicz, Lasota & Igumenshchev 2000). If outflows are as 
strong as Blandford & Begelman (1999) suggest, both jet mecha- 
nisms we have discussed here are impossible. One mechanism de- 
pends directly on Madaf > while the other depends on the central 
pressure ptot which is roughly proportional to Madaf- If there are 
heavy outflows, there is just not enough energy to power the ob- 
served jets. 
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Figure 5 shows the variation of SS with r for the solutions 
shown in Figs. 1 and 3. We have normalized by the local gravi- 
tational potential to obtain a dimensionless measure of how 
unbound the gas is at each radius. All the solutions start with a large 
value of ,i^/|<5| outside the Bondi radius, but this just means that 
the gas out there has a finite thermal energy but very little gravita- 
tional binding energy. As the gas flows in, .'M hardly changes while 
$ increases rapidly in magnitude, so the ratio S§ /\^\ decreases 
rapidly. In the case of the Bondi-like solution with ^ = 0.11, the 
decrease continues all the way down to the horizon and there is 
no tendency to form an outflow. The ADAF solutions, on the other 
hand, follow the Bondi solution down to r ~ 10^''^, after which 
viscous dissipation causes ,^/|<I>| to increase up to a maximum 
value of about 0.1. We do not know if this value is large enough 
to strongly suppress the mass accretion rate on to the black hole. 
However, since the Bernoulli parameter is large only over a limited 
range of radius, it is conceivable that outflows reduce the mass ac- 
cretion rate in the ADAF solutions by no more than a factor of a 
few, rather than by orders of magntiude. In this case, both of our jet 
mechanisms are likely to survive. 

Our model assumes a single point source of gravitation, 
whereas a realistic galactic nucleus has significant stellar mass from 
the inner galaxy and any nuclear star cluster This will make the 
outer gas more gravitationallly bound than appears from Fig. 5. 
Radiative cooling (which is ignored in our model) can also have an 
effect. Using the gas properties at the Bondi radius given in Allen et 
al. (2006) for M87 we find that the radiative cooling time is about 
2.5 percent of the flow time (r/cs). The ratio of these timescales 
increases at smaller radii, making radiative cooling unimportant 
there. Compton cooling however increases with decreasing radius 
and may be important at precisely those inner radii where SS shows 
an increase in Fig. 5. Thus, the importance of the Bernoulli param- 
eter may be reduced even further 

Another important effect pointed out by Narayan & Yi (1994) 
is that viscous dissipation, coupled with the lack of radiative cool- 
ing, causes the entropy of the gas in an ADAF to increase in- 
ward, making the flow convectively unstable by the Schwarzschild 
criterion. Convective effects have not been included in the one- 
dimensional model we have considered in this paper. Narayan, 
Igumenshchev & Abramowicz (2000) and Quataert & Gruzinov 
(2000) discussed the physics of convection-dominated accretion 
flows (CDAFs) and concluded that such flows would differ enor- 
mously from ADAFs. In particular, if one considers a self-similar 
model, the density, pressure and mass accretion rate of a CDAF, as 
measured at the black hole, are predicted to be a factor ~ rg/rg ~ 
10"^ (for our models) times the corresponding values for an ADAF 
with the same outer boundary conditions. Even if the real effect is 
only a small fraction of this analytical prediction, it would reduce 
jet power to a level far below that observed. 

Numerical hydrodynamic simulations confirm the presence 
of convection in ADAFs (Stone et al. 1999; Igumenshchev & 
Abramowicz 1999, 2000; Igumenshchev, Abramowicz & Narayan 
2000) and indicate substantial suppression of the mass accretion 
rate into the black hole. This is problematic for jet production. 
However, the relative importance of ADAFs vs CDAFs has been 
debated (e.g., Igumenshchev & Abramowicz 1999; Abramowicz et 
al. 2002; Lu, Li & Gu 2004), and the issue is still unresolved. 

A measure of convective instability is the Brunt- Vaisala fre- 
quency A' which for a spherically symmetric system is given by 

p dr dr 




0246 3 46 

log(r/rg) log(r/r|,) 



Figure 6. Left: Variation with radius of the quantity —N-t1^^, which is a 
measure of convective instability, for two of the four solutions shown in 
Fig. 1. (The other two are below the bottom of the plot.) The dotted line 
indicates our nominal threshold for becoming convection-dominated. The 
more rapidly spinning solution (^ = 85) goes above the dotted line over a 
range of radius, and we expect it to become a CDAF in this zone. The less 
rapidly rotating solution (^ =12) lies below the dotted line and probably 
will not become a CDAF. Right: Corresponding results for the solutions 
shown in Fig. 3. The two solutions with a = 0.316 and 0.1 lie below the 
dotted Hne and will not be convection-dominated, while the solutions with 
smaller values of a are expected to become CDAFs. Note that none of the 
curves in the two panels extends inside the sonic radius. This is because 
we have assumed the flow to be inviscid once it becomes supersonic. The 
specific entropy is then constant ( ^2.2) and = 0. 

A system is convectively unstable if A^^ < 0. When this happens, 
the quantity |A'| = V—N^ measures the growth rate of the insta- 
bility. We imagine that convection is important and takes over the 
dynamics of the flow only when the growth time scale of the in- 
stability is shorter than the accretion time /acc = whereas in 
the opposite limit, we expect convection to be a minor perturba- 
tion. Motivated by this argument, we show in Fig. 6 the profiles of 
the dimensionless quantity —N^t^^c versus r for a selection of our 
numerical solutions. Most of the solutions of interest to us, viz., 
those with relatively large values of a and small rotation parame- 
ters have —N'^tlcc < 1- Convection is probably unimportant in 
these cases. The results shown in the panel on the right agree with 
the a trend discussed by Igumenshchev & Abramowicz (1999) and 
Lu et al. (2004). 

In the discussion so far, we have ignored the effect of mag- 
netic fields. Magnetohydrodynamic (MHD) simulations of ADAFs 
(Stone & Pringle 2001; Hawley, Balbus & Stone 2001; Igumen- 
shchev & Narayan 2002; Igumenshchev, Narayan & Abramow- 
icz 2003; Pen, Matzner & Wong 2003; Igumenshchev 2006, 2008; 
Pang et al. 2010) have not so far exhibited strong unbound winds, 
but they nevertheless have shallow density profiles and reduced 
M. The flows exhibit vigorous turbulence and they transport en- 
ergy outward, just as one expects with a convective flow. However, 
whether or not the turbulence can be described as convection is 
unclear (Narayan et al. 2002, Balbus & Hawley 2002; Pen et al. 
2003). 

From the point of view of our present study, the relevant ques- 
tion is how much is the mass accretion rate suppressed relative to 
the Bondi rate as a result of MHD effects. The best scalings from 
current simulations suggest that the accretion rate, and therefore the 
jet power, is reduced by at least three orders of magnitude. More- 
over, M is found to be reduced substantially even in the case of a 
non-rotating Bondi flow (Igumenshchev & Narayan 2002; Igumen- 
shchev 2006). This last result is particularly worrisome in view of 
the results shown in Fig. 6. There it appeared that, so long as a 
is relatively large and the rotation is small, hydrodynamic convec- 
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tion and the consequent suppression of M are not an issue. How- 
ever, it appears that magnetic fields completely alter the situation 
and strongly suppress even Bondi accretion. This effect needs to be 
confirmed with further studies and the physics of the phenomenon 
needs to be identified. 

Apart from convection, thermal conduction might also mod- 
ify the dynamics of radiatively inefficient accretion (e.g., Johnson 
& Quataert 2007; Sharma, Quataert & Stone 2008; Shcherbakov & 
Baganoff 2010). Numerical MHD simulations including conduc- 
tion and covering an adequate range of radius are yet to be carried 
out. 

If any of the effects discussed here succeeds in strongly sup- 
pressing the mass accretion rate in quasi-spherical accretion flows, 
we would be left with the puzzle of why the observed jet power Pj 
in many nearby low-luminosity galactic nuclei tracks the estimated 
Bondi mass accretion rate Mb i^H'- Vj = Pj/M^c^ ~ 2%. The ob- 
served jets would require a different explanation that is not related 
to a hot accretion flow, and the apparent correlation with properties 
of hot gas at the Bondi radius must be a coincidence. 

A disk of dense cool gas emitting optical lines is often seen in 
the nuclei of elliptical galaxies. Macchetto et al. (1997) have stud- 
ied the gas distribution in M87 and used the kinematics of the gas 
to estimate the mass of the black hole. The result is slightly less 
than that now determined from stellar kinematics by Gebhardt & 
Thomas (2009). The gas disk has a central hole a few pc in size. 
Line widths are large, possibly indicating the action of nongravita- 
tional forces in this and other elliptical galaxies (Verdoes-Kleijn, 
van der Marel & Noel-Storr 2006). We do not think that small 
masses of such cool gas are incompatible with the existence of a 
giant ADAF. 

Finally, we note that estimates showing that Bondi accretion 
cannot yield a sufficient mass accretion rate (e.g., Rafferty et al. 
2006) are based on the temperature and density values inferred at 
radii far outside r^. Small quantities of cooler gas at the centre, say 
at 0.7 keV instead of 3 keV, which is in pressure equilibrium with 
the surrounding gas would have a correspondingly higher density 
and allow a much higher Bondi flow rate (Mb °= T^^l'^). Such cool 
gas is often observed in nearby clusters where the innermost re- 
gions are spatially-resolved. The rate will in practice be even higher 
than this simple scaling as the pressure will be higher at the Bondi 
radius due to the weight of intervening gas. 



5 SUMMARY 

We have shown in this paper that accretion can occur from a hot 
atmosphere at close to the Bondi rate, provided the external gas 
near the Bondi radius rotates relatively slowly (less than a few tens 
of percent of the Keplerian rate) and the viscosity parameter a is 
fairly large (> 0.1). The non-radiative numerical ADAF solutions 
computed here may be relevant to the nuclei of massive elliptical 
galaxies hosting a central black hole surrounded by a hot gaseous 
atmosphere. The mass accretion rate is large enough that it requires 
only a small fraction (~ 2 — 5%) of the accretion energy to power 
the observed jets in nearby elliptical galaxies. At the same time, the 
jets can heat the surrounding gas and prevent the hot atmosphere 
from radiatively cooling and collapsing into the centre. This feed- 
back mechanism could also have an effect on the evolution of the 
galaxy itself. These results require that mass outflows (other than 
that associated with the jets), convective energy transport and MHD 
effects are weaker than usually assumed. 
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